Arbitrarily Accurate Dynamical Control in Open Quantum Systems 
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We show that open-loop dynamical control techniques may be used to synthesize unitary trans- 
formations in open quantum systems in such a way that decoherence is perturbatively compensated 
for to a desired (in principle arbitrarily high) level of accuracy, which depends only on the strength 
of the relevant errors and the achievable rate of control modulation. Our constructive and fully ana- 
lytical solution employs concatenated dynamically corrected gates, and is applicable independently 
of detailed knowledge of the system-environment interactions and environment dynamics. Explicit 
implications for boosting quantum gate fidelities in realistic scenarios are addressed. 
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The demand for an exquisite degree of control over 
the dynamics of open quantum systems is widespread 
across quantum physics and engineering, ranging from 
high-resolution spectroscopy and chemical reaction con- 
trol [1], to quantum- limited metrology [2] and quantum 
information processing (QIP) Achieving a sufficiently 
small 'error per gate' (EPG) is, in particular, an essen- 
tial ingredient to ensuring that fault-tolerant quantum 
computation is possible in spite of the decoherence that 
inevitably plagues real-world devices. While closed-loop 
techniques, in the form of fault-tolerant quantum error 
correction (QEC) [4, 5], offer thus far the only complete 
prescription to meet this challenge, open-loop dynamical 
QEC is emerging as a promising alternative. Inspired 
by coherent averaging in magnetic resonance [G] and ex- 
emplified in its simplest form by dynamical decoupling 
(DD) [7"')], dynamical QEC aims to suppress the interac- 
tion between the system and its quantitatively unspecified 
environment through suitable sequences of unitary oper- 
ations. Recently, DD has enabled dccohercncc-protectcd 
storage in QIP platforms as diverse as electron-nuclear 
systems [lU], photonics qubits [11], and trapped ions [12], 
as well as found application in suppressing collisional de- 
coherence in cold atoms [13]. 

As the gap between theory and implementations 
shrinks, and a growing experimental effort is devoted 
to robust manipulation of quantum states, it is imper- 
ative that realistic constraints be accommodated from 
the outset in dynamical QEC design. In practice, con- 
trol resources always entail finite power and bandwidth, 
thus precluding instantaneous ('bang-bang' [7]) pulses. 
A path toward dccohercncc-protectcd unitary operations 
was recently proposed based on dynamically corrected 
gates (DCGs) [14]. A DCG may be viewed as a com- 
posite quantum gate constructed from individual ('prim- 
itive') building blocks whose errors combine non-linearly 
to achieve a substantially smaller net error. If Tmin is 
the minimum duration over which each primitive gate is 
effected ('switching time') and \\Hf.\\ the strength of the 
error-inducing Hamiltonian, DCGs remove the effect of 



Hf. to the leading (first) order, that is, the resulting EPG 
scales as ©[(rminll-f^ell)^]- This prompts the following key 
question: Can one make DCGs as accurate as desired, 
using realistic control resources? The answer is not obvi- 
ous. Schemes capable of arbitrarily suppressing decoher- 
ence during storage have been identified - notably, con- 
catenated DD [n] and recent optimized protocols [9, 15] - 
but, thus far, only in the bang-bang limit. With bounded 
controls, decoherence suppression up to the second or- 
der (with leading corrections 0[(ri„in||i/e|j)'^]) may be 
achieved by using a time-symmetrized Euler DD (EDD) 
protocol [IG] (see also [17]) - however, this procedure ex- 
tends neither to generic quantum gates nor to generic 
open quantum systems [14, 18]. 

Here we show that decoherence suppression can in 
principle be pushed to an order limited only by the 
strength of the relevant errors and the achievable rate 
of control modulation. We do this by combining DCG 
constructions with recursive design - resulting in con- 
catenated DCGs. While perturbative in nature, our so- 
lution is fully analytical, laying the foundation for rigor- 
ous complexity analysis and optimization in dynamical 
QEC. Not only do concatenated DCGs exist for arbi- 
trary finite-dimensional open quantum systems with a 
bounded U-ffell, but they are also highly portable, in the 
sense that no quantitative knowledge of the underlying 
interaction Hamiltonian is assumed, beyond its algebraic 
form. Since arbitrarily accurate open-loop compensation 
techniques for classical (static) control errors are known 
[19], our results imply that no fundamental limitation 
arises due to a quantum (dynamic) bath. From a practi- 
cal standpoint, concatenated DCGs offer the first system- 
atic feedback-free framework for designing quantum gates 
which can achieve the arbitrarily high levels of protection 
against decoherence demanded by high-fidelity quantum 
control, and in particular QIP. 

Control-theoretic setting. — Let S be the target quan- 
tum system, coupled to its quantum bath B via an inter- 
action Hamiltonian Hsb, with respective Hilbert space 
Tis and Hb- Wc assume that the total error Hamilto- 
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nian may be described by a joint time-independent 
operator of the following form: 

a 

where the contribution Hs^e accounts for undesired 'drift' 
terms possibly present in the system's internal Hamilto- 
nian, {Sa} is a Hermitian basis of operators acting on 
Hs, and Ba are bounded (potentially unknown) oper- 
ators acting on Hb- The vector space fig spanned by 
{Sa ® Ba}, with non-zero Sq,, defines the error model 
and is uniquely determined by the system components 
{5'q,} appearing in the expansion of in Eq. (1). In 
dynamical QEC, a classical controller is adjoined to S 
through a time-dependent Hamiltonian Hctriit). Several 
constraints may restrict, in reality, the degree of con- 
trol that is available through Hctiiit). In particular, we 
account for realistic control modulations to be bounded 
in amplitude and spectral bandwidth by requiring that 
ll-ffctri(i) II < oo and Tmin > 0, respectively. In an ideal 
situation where Hsb = = -ffs^e, the errors arise from 
imperfections in the controller only. Our goal here is to 
present a proof-of-concept for robustness against decoher- 
ence errors during gates, thus we ignore such contribu- 
tions henceforth. We furthermore assume that universal 
control over S is achievable [3] (see [14] for departures 
from this 'minimal' setting). 

Consider a unitary operation (gate) Q on S, which is 
ideally realized by letting Hcti\{t) — HQ{t) over an inter- 
val [ti,ti + t]. We use Q[t] or Q to denote this imple- 
mentation of Q when there is no ambiguity. If H^. ^ 0, 
the joint propagator is J7q[t] = Texp[—i J^^'^'^ {HQ(t) + 
He)dt] {h = 1 and T denotes time-ordering). The effect 
of He may be isolated through a Hermitian 'error action 
operator' i?Q[T] [14], where Uq^^-^ = Q exp(— zE'q [,.]). Let 
the initial states of S and i? be = \4'){'4'\ and pB- Ide- 
ally, the action of Q would result in P%{t) = Qp%Q\ 
while B evolves independently. In contrast, the ac- 
tual evolution is coupled and ps{t) = TrB{UQ[T]P% ® 
PbUq^t]}- appropriate performance measure for con- 
trol is the trace-norm distance A = WpsiT) — P^s^''')\\'^ ['"']■ 
One can show that A < ||modB(ii'Q[^])||, independently of 
the initial state, where modB(i?) = E—t^^^^Is®^t:s{E) 
is a projector that removes the pure-bath terms in E 
[20, 21]. Thus, rj = ||modB(i?Q[T-])|| may be taken to 
quantify the resulting EPG. 

Given a desired unitary gate Q, our task is to syn- 
thesize a control modulation that approximates Q with 
an EPG scaling as 0{t^^^), for an arbitrary positive 
integer ^ and compatible with the stated control con- 
straints. Any such construction defines an ^-th order 
DCG, DCG^^^ "Naive" implementations that simply cor- 
respond to turning on, say, a constant Hctrxit) for a dura- 
tion T yield an EPG that scales (approximately) linearly 
with T in the presence of He and arc included as £ = 0. 



Primitive gates are correspondingly denoted by {Q'*^'}. 
Free evolution under He [i?ctri {t) = 0] may be viewed also 
as a zeroth order 'no-operation' (noop) gate. Consider 
a combined gate {Qi ■ ■ • Q n)^i^ ^i] formed as a sequence 

of N primitive gates [tjv] • • • (52°' ['''2]Qi'' ['''i] applied 
back to back. We may compute the total error for the 
combined gate using 

■^(Qi---Qiv)Er.] = ^Qi[ri] + (2) 

where Pj = Qj ■ ■ - Qi, j = 1, . . . , N — 1, is the 'partial' 
control propagator at the end of the jth segment [22] and 
\\E2+\\ = ||[^Q.[T.]:^Qj[rj]||) includes the leading 

second-order corrections. Note that individual EPGs do 
not simply add up, but are 'modulated ' by the control tra- 
jectory: in dynamical QEC, this provides the basic mech- 
anism enabling the effect of He to be pcrturbativcly can- 
celed, without quantitative knowledge of the bath. 

Concatenated DCGs: Construction and performance 
bound. — Our strategy to achieve an arbitrary order of 
cancellation is to invoke a recursive construction, that 
generates {£ + l)"^-order gates using £"^-order building 
blocks, close in spirit to concatenated DD [S]. In what 
follows, given a target gate Q, we shall term an i-th or- 
der balance pair as any pair of gates ( Jg' , (5* ' ) whose 
errors coincide up to the leading (£ -I- l)**^-order, that is, 
mods {Ej^ ) = mods {Eq,) + 0(r;^+f ) . ^ 

The first step in the recursion requires demonstrat- 
ing that, given Q, a corresponding DCG^^l can be con- 
structed out of the available universal set of primi- 
tive gates. As established in [14], this necessitates two 
main ingredients: (i) An EDD protocol for generating a 
DCG^^l-implementation of NOOP; (ii) A zeroth-order bal- 
ance pair {I^Q , Q'^*^), with the same leading order error 
Eq, belonging to a 'correctable error space' Q,ti 3 fig. In 
order to achieve (i), a set of unitary gates V = {Di}f^^ 
is identified, such that the map nx)[i?] = ^ X^iLi Di^Di 
'decouples' all errors E in Hu: modB(nD[i?]) = [7]. In 
EDD [Ki], T> represents a group Q (faithfully and pro- 
jectively), with order \Q\ = d and a set of genera- 
tors r ~ {Fj}. Let G{Q,T) denote the Cayley graph of 
Q with respect to F. The required EDD sequence (of 
length Ni = dm) is constructed by consecutively apply- 
ing the generators Fj as gates, in the order determined 
by an Eulerian cycle on G{G,^), starting (and ending) 
at the identity vertex. Is- Suppose that each Fj is im- 
plemented as a zeroth order gate with Ep. G Vlx). Then, 
using Eq. (2) and the decoupling property of G, the net 
error 77EDD = 0{ma.Xj WEp^W^) = C(T^i„), as desired. 

In order to extend the EDD construction to a gate Q 
other than NOOP, some information about how primitive 
gates are implemented is required, due to a NoGo theo- 
rem for 'control-oblivious' design proved in [14]; balance 
pairs are required as ingredient (ii) precisely for this pur- 
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pose. Once such a balance pair {Iq , Q* ) is found, a 
DCG^^l sequence for Q is obtained by 'augmenting' the 
Cayley graph for noop: /g' gates arc inserted in the 
EDD sequence at points where the corresponding Eule- 
rian path visits each non-identity vertex for the first time, 
and finally Qp' is applied after the last Is vertex. The 
resulting error for the combined sequence is then given by 
||modi3(-BQ[ii)|| = WiRodBiEEBB +'nviEQ,))\\ = plus 

corrections given by error — 0{t^^-^^), as desired. 

Our balance pair construction "strikes a balance" by 
stretching control profiles. Not only does this ensure a 
fully portable recipe in terms of control inputs (note that 
previous DCG constructions [14], required access to sign- 
reversed control profiles), but it is in fact crucial for pro- 
ducing balance pairs of high-order gates. In practice, 
a stretched control input results in slower gates. More 
concretely, for a gate stretched by a factor r, the gat- 
ing Hamiltonian is "stretched" by a factor r: HQ{t) (-t- 
HQ\ti + r{t — ti)]/r; thus the same target gate Q is ap- 
proximated, but with a different EPG. The stretched gate 
is symbolically indicated by Q[rT] but all gate construc- 
tions are still subject to the original bandwidth/power 
constraints. The general recipe for balance pairs of any 
order £ is given in terms of Q^^l and (^"^ order 

implementations of Q and its inverse) by the following 
pair of composite gates [21]: 

jM = Q-i.M[r]QM[2V(^+i)r], 
Ql''[r] = QM[r]g"i-M[r]QM[r]. (3) 

Note that i'q implements the NOOP gate (Is) over a du- 
ration of T + 2^/(^"''^V. while Qi^' implements Q over 3r. 
Setting £ = thus completes the construction of f^'-order 
gates using 0*^-order building blocks. 

A (universal) set of gates {Q[^1[t£]} with EPG = 
0{t^^) can be constructed recursively for arbitrary 

£ > 1 at this point. Let fig denote the error model for 
all £*'^-order gates: EQ[t] S ri[f' [23]. We can then iden- 
tify the smallest DD group Q^^^ that decouples fl^eK with 
m£ group generators r^^l ~ {F^^^}, and a correspond- 
ing Cayley graph G{Q^^\T^^^). We modify this graph 
by attaching self-directed edges representing Iq to all 
vertexes except /g, and add a new vertex Q by con- 
necting it to the /^-vertex through an edge representing 
Qi . By construction, every edge in this graph repre- 
sents a DCGM. We now proceed as in first-order DCGs, 
and implement the sequence QI^^^I [Ti+i] by following the 
Ni = dimi + df edges of an Eulcrian path on the mod- 
ified graph for G(a[^I,rM), starting at Is and stopping 
at Q, and applying the corresponding Q^^ gates back to 
back. If the latter are implemented with duration (be- 
fore stretching), then the combined total duration r^+i 



satisfies: 

Te+i = [dime + {dt - 1)(1 + 2^''^^+^'^) + 3]t,. 

By iterating, we obtain a cascade of stretched primitive 
control profiles of duration T( < {xeY'''o^ where = 
di{mi -\- 3) [24]. Starting with the primitive gates [£ = 0] 
of duration tq, the above construction generates a DCG^^l 
for any Q and £, with a net error upper bounded by: 

^^DCG < c(X£)'Vo \\HsB+Hs.e\\ (4X^^0 ||i^e!|)' , (4) 

where c~0{\). While a detailed proof will be presented 
elsewhere [21], two main steps are involved. First, we 
prove that for any gate Q[^'[t], owing to the recursive 
design of the sequence, the Magnus expansion of the error 
Eqye] in terms of the toggling frame error Hamiltonian 

He{t) = Ul^.^^{t)IIJJctr\{t) contains only terms that start 
at 0(t^+1) (modulo the pure bath terms). Second, we 
bound these higher order terms (hence ]|modB(i?Q[«] ]|)) 
using standard operator inequalities. 

Concatenated DCGs: Analysis and applications. — 
From a practical perspective, our result above shows 
that concatenated DCGs offer concrete error reduc- 
tion over primitive gates and DCG^^l: for any fixed 
achievable switching time tq, minimizing the bound in 
Eq. (4) yields an optimal concatenation level io-pt = 
[— ^(log^(4)]/Je|)To) + 1)J, where [.tJ is the largest in- 
teger < a;, and x = Xfopt ■ Substituting ^opt into Eq. (4) 
yields an EPG bound that may enable scalable QIP even 
if primitive gates, whose EPG is given by \\HsB+Hs,e\\To, 
are above the accuracy threshold of (non-Markovian) 
fault-tolerant QEC [25]. Ultimately, the viability of a 
concatenated DCG will be dictated by system-dependent 
implementation trade-offs, based on both the total gate 
duration and minimum switching time. 

To illustrate our general construction, consider the 
paradigmatic case of a single qubit undergoing arbitrary 
decoherence, whereby filf' = O'^^ = spanjcTQ (Si Ba}, and 
(Tq = {/, X, Y, Z} are the identity and the Pauli matri- 
ces. The corresponding DD group 5 = Z2 x Z2 is (projec- 
tively) represented as {I, X,Y, Z}, and is generated by, 
e.g., r = {X,Y}. EDD is given by XYXYYXYX [10]. 
Starting with a Q[^l[r£] gate, we have: 

in which /g' and q\^'^ arc defined in Eq. (3), and are 
themselves combinations of stretched Q^^' and Q~^'^^^. 
The length of the combined new sequence is given by 
T£+i = [14 + 3 X 2i/(^+i)]r£ < 20t£ = xn- Representa- 
tive simulation results are presented in Fig. 1. Param- 
eters have been chosen to mimic a high-quality silicon 
(Si) quantum dot, where the electron spin qubit under- 
goes hyperfine-induced decoherence due to a fraction ?» 1 
ppm of non-zero spin ^^Si nuclei (about one order of mag- 
nitude larger than currently achieved isotopically purified 
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FIG. 1: Fidelity of increasing DCG orders for Q = 
exp[-i(7r/3)X] applied to = (|0) + |l})/ %/2, and ps fully 
mixed. The fidelity is / = Try^^/pi^pt^/pZ, where pt (pa) is 
the target (actual) final qubit state, and satisfies 1 — A < / < 
y/1 — [28]. The bath spins couple to the central system 
spin S via a Heisenberg interaction Hsb = X^i^iJ*^ ' > 
with ji randomly picked in [0, J], J = 10 (arbitrary units). 
The bath spins evolve under a dipolar interaction Hb ~ 
bi]{Ix'' ^x^ +Iy -'^^z^ ^z^)^ with bij randomly picked 
in [0, 10~^]. Primitive gates were implemented using rectan- 
gular pulse shapes to allow for numerically exact simulations. 
For fixed J, reducing Tmin may be understood in terms of a 
finer temporal resolution of the resulting 'digitized pulse pro- 
file'. Note that at logio(Tmin J) = -5.5 the fidelity for DCG'^' 
is 13 orders of magnitude better than for DCG'"', in spite of 
the gate taking nearly 6500 times longer. 

Si [2(i]). Consistent with the perturbative nature of the 
error cancellation in dynamical QEC, the improvement 
due to the recursive design is manifested in the increasing 
slopes associated with higher concatenation levels once 
the gating time Tmin is sufficiently short. Concatenated 
DCGs may also prove instrumental in reducing gating 
errors on a recently proposed logical qubit encoded in 
the singlet/triplet spin manifold of a Si double quantum 
dot [27], in particular to protect exchange-based logical 
Z rotations against magnetic field noise. Finally, we ex- 
pect the same control sequences to be effective for non- 
Markovian decoherence induced by unbounded environ- 
ments with a sufficiently 'hard' spectral cutoff [9]. 

Conclusion. — We have presented a general construc- 
tive solution to the problem of generating arbitrarily ac- 
curate quantum gates with finite control resources in an 
open-loop setting. In addition to settling a fundamental 
question, our results point to several venues for further 
investigation. On the theory side, improved construc- 
tions should incorporate more realistic local baths and 
combinatorial Eulerian design [29]. From an implemen- 
tation perspective, making contact with optimal-control 
formulations [30] may allow to boost efficiency in exper- 
imentally available control platforms. 
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